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Abstract 

Diffusion  of  oxygen  through  the  thin  liquid  film  (TLF)  covering  the  surface  of  the  negative  active  mass  (NAM)  crystals  and  impeded  charge 
transfer  through  the  electric  double  layer  (EDL)  were  assumed  to  be  the  rate  limiting  stages  of  the  closed  oxygen  cycle  in  VRLAB.  Based  on 
the  theory  of  colloid  systems,  the  thickness  of  the  TLF  was  determined  and  its  dependence  on  pore  radius  and  surface  tension  at  the  gas/ 
electrolyte  interface  was  estimated.  The  initial  and  boundary  conditions  were  defined  and  Fick’s  equation  was  solved  for  this  case.  The 
diffusion  kinetics  results  obtained  explain  the  dependence  of  the  recombination  current  on  the  metal  surface  on  which  the  electrochemical 
reaction  takes  place,  the  thickness  of  the  TLF,  the  concentration  of  oxygen  at  both  film  surfaces,  and  the  time  of  polarization.  Using  the  film 
thickness  equation  derived  by  the  colloid  chemistry  theory,  a  general  equation  was  developed  which  determines  the  oxygen  recombination 
rate  dependence  on  the  thin  liquid  film  parameters,  NAM  structure  and  polarization  time.  The  validity  of  this  equation  was  experimentally 
proved  for  low  and  high  electrolyte  saturation  levels  of  the  plates,  as  well  as  for  negative  plates  produced  employing  various  technologies.  This 
equation  makes  it  possible  to  draw  conclusions  about  the  possible  ways  of  improvement  of  the  efficiency  of  the  oxygen  cycle  and  about  the 
technologies  for  negative  plate  production  for  VRLAB  applications. 
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1.  Introduction 

The  closed  oxygen  cycle  (COC)  comprises  the  following 

elementary  processes: 

(a)  Oxygen  evolution  on  the  surface  of  the  positive  plate. 

(b)  Transport  of  02  through  the  AGM  separator  to  the 
negative  plate.  Oxygen  transport  proceeds  through  gas 
channels  in  AGM  when  the  electrolyte  saturation  in  it 
is  <96%. 

(c)  The  lead  particles  in  the  negative  plate  are  covered  by  a 
layer  of  electrolyte.  The  gas  (02)  must  dissolve  in  this 
layer  and  the  oxygen  molecules  must  diffuse  to  the  lead 
surface  at  which  the  reduction  proceeds  through  accept¬ 
ing  electrons.  The  concentration  of  oxygen  in  the 
sulfuric  acid  solution  is  about  8  x  10-7  mol  cm-3.  The 
diffusion  coefficient  of  02  is  about  1— 2  x  10-5  cm  s-1 
(depending  on  acid  concentration).  This  figure  shows 
that  the  liquid  layer  covering  the  metal  surface  may  slow 
down  markedly  the  rate  of  the  COC. 
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(d)  Oxygen  may  react  at  the  electrode  surface  following 


two  different  mechanisms: 
o  Chemical  one 

Pb  +  j02  — >  PbO  (a) 

PbO  +  H2S04  PbS04  +  H20  (b) 

PbS04  +  2H+  +  2e“  ->  Pb  +  H2S04  (c) 

o  Electrochemical  one 

z02  +  4H+  +  4e“  -►  2H20  (d) 


(e)  H20  transport  from  the  negative  towards  the  positive 
plate. 

The  slowest  elementary  process  of  the  COC  has  been  the 
object  of  several  studies: 

•  It  has  been  assumed  that  the  slowest  stage  of  the  COC  is 
the  diffusion  through  a  liquid  layer  wetting  the  surface  of 
the  negative  plate  [1-4].  The  authors  investigate  the  case 
of  constant  gas  volume  and  cell  temperature.  The  pressure 
in  the  cell  will  increase  until  the  valve  opens.  If  the 
charging  current  is  sufficiently  high  the  gas  pressure  in 
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the  cell  in  this  case  will  be  kept  constant  by  the  valve.  On 
long  overcharge  when  the  state  of  charge  (SOC)  is  100%, 
the  composition  of  the  outgoing  gas  will  rapidly  reach  a 
H2:02  ratio  of  2:1  [4]. 

•  It  has  been  established  that  the  recombination  rate 
depends  but  very  slightly  on  the  potential  of  the  negative 
plate.  The  process  of  02  reduction  at  rotating  disk  elec¬ 
trodes  is  a  diffusion  limited  one,  it  can  be  represented  by 
the  Levich’s  equation  [5-8]; 

•  The  influence  of  H2S04  concentration  on  the  kinetics  of 
02  reduction  on  lead  electrodes  was  determined  in  [9]; 

•  The  transport  of  02  to  the  surface  of  the  negative  plate  and 
the  02  reduction  process  were  investigated  by  mathema¬ 
tical  models  and  computer  simulations  proposed  by  Ber- 
nardi  and  Carpenter  [10],  Newman  and  Tiedeman  [11],  and 
Gu  et  al.  [12].  An  important  assumption  in  these  computer 
simulations  is  the  consideration  of  impeded  charge  transfer 
through  the  electric  double  layer  (EDL)  during  oxygen 
reduction  on  the  negative  plate.  This  assumption  improves 
markedly  the  match  between  the  calculated  and  the  experi¬ 
mental  data  in  all  three  investigations  [10-12].  The  authors 
use  different  variants  of  the  general  Volmer’s  equation  for 
describing  the  kinetics  of  oxygen  reduction. 

•  A  mechanism  of  the  processes  taking  place  during  oxygen 
evolution  on  Pb02  was  proposed  in  [13,14]. 

The  properties  of  the  thin  liquid  film  (TLF)  of  electrolyte 
covering  the  Pb  surface  and  their  impact  on  oxygen  diffusion 
were  not  considered  in  the  above  studies.  The  aim  of  the 
present  investigation  is  to  determine  the  influence  of  the 
properties  of  the  thin  liquid  film  on  the  kinetics  of  oxygen 
diffusion  and  its  reduction  on  the  Pb  surface. 


2.  Oxygen  reduction  processes  taking  place  on 
the  lead  surface 

We  took  the  following  features  into  account  when  con¬ 
sidering  the  oxygen  diffusion  and  reduction  processes  at  the 
Pb  surface: 

(a)  The  negative  plate  is  characterized  by  high  porosity  and 
large  surface  area  (0.8  m2  g-1).  Under  VRLAB  operat¬ 
ing  conditions,  the  electrolyte  saturation  of  the  negative 
plate  is  less  than  90%.  The  gas  (02)  will  penetrate  into 
the  volume  of  NAM  through  the  open  gas  macro  pores 
of  the  plate. 

(b)  The  thickness  of  the  thin  liquid  film  (TLF)  depends  on 
the  capillary  pressure  in  the  open  gas-liquid  pores. 
Under  the  action  of  the  capillary  pressure  the  film 
covering  the  surface  of  the  lead  crystals  in  the  PAM 
becomes  thinner. 

(c)  Water  is  formed  as  a  result  of  oxygen  reduction.  Thus, 
there  is  a  tendency  for  the  thickness  of  the  TLF  to 
increase  depending  on  the  rate  of  02  reduction. 

(d)  The  diffusion  transport  of  oxygen  follows  Fick’s 
equation  and  depends  on  the  diffusion  coefficient  of 


oxygen  in  the  H2S04  solution,  on  the  concentration  of 
oxygen  on  both  phase  boundaries  of  the  TLF  and  on  the 
thickness  of  the  TLF. 

(e)  Oxygen  diffusion  through  the  TLF  and  charge  transfer 
through  the  electric  double  layer  (EDL)  are  the  slowest 
steps  of  the  process  of  reduction  of  oxygen  molecules. 

(f)  NAM  contains  expanders  (polymeric  surfactants) 
which  have  an  influence  on  the  TLF  and  especially 
on  the  electric  double  layer  at  the  metal  surface.  Thus 
expanders  can  influence  the  processes  of  oxygen 
reduction. 


3.  Evaluation  of  the  TLF  thickness  and  its  influence 
on  COC  efficiency 

The  thickness  of  the  thin  liquid  film  h  is  a  dynamic 
parameter  which  depends  on  both  the  changes  in  capillary 
pressure  in  the  pores  (tending  to  reduce  h)  and  on  the 
increased  rate  of  recombination  (generation  of  water,  tend¬ 
ing  to  increase  h).  The  TLF  thickness  depends  also  on  the 
nature  of  the  expanders  used  which  determine  the  surface 
tension,  a  on  the  phase  boundary  gas-liquid.  Under  the 
influence  of  the  above  parameters  the  thickness  of  the  TLF 
will  tend  to  reach  a  certain  steady  state  value.  The  latter  will 
depend  also  on  the  electrolyte  saturation  level  in  the  negative 
plate.  The  thickness  can  be  estimated  using  the  theoretical 
methods  of  colloid  chemistry. 

An  open  gas  macropore  in  the  negative  plate  is  diagram- 
matically  represented  in  Fig.  1.  Such  macropores  are  formed 
as  a  result  of  reduced  electrolyte  saturation  in  the  negative 
plate.  The  capillary  properties  of  the  water  solution  open  the 
pore  and  make  its  cross-section  round  shaped.  Let  us  assume 
that  a  wetting  film  of  electrolyte  is  formed  between  the  pore 
wall  and  the  lead  crystal  surface  (Fig.  2). 

A  wetting  film  with  thickness  h  is  formed  in  a  pore  with 
radius  R.  The  balance  of  forces  on  the  curved  meniscus  area 
determines  the  capillary  pressure  Pa  in  the  pore,  which  is 
equal  to  PCT  =  P\  —  Pq.  P i  is  the  pressure  in  the  gas  phase 
and  P0  the  pressure  in  the  liquid  phase.  The  force  balance  on 
the  gas/liquid  interface  of  the  film  determines  the  disjoining 
pressure  77  :  77  =  Pi  —  P0.  At  equilibrium  =  77. 

The  stability  conditions  for  such  a  wetting  film  can  be 
calculated  from  the  viewpoint  of  the  so  called  DLVO  theory 
( Derjaguin ,  Landau ,  Verwey  and  Overbeek )  [15-19].  The 
stability  of  the  film  depends  on  the  sign  of  the  disjoining 
pressure  77,  which  is  a  basic  thermodynamic  characteristic 
of  the  TLF.  77  accounts  for  the  force  of  interaction  between 
the  two  surfaces  of  the  film. 

The  DLVO  theory  considers  77  as  a  superposition  of  two 
basic  universal  molecular  interactions,  the  electrostatic  and 
the  van  der  Waals  ones: 

77  =  77ei  +  77vdw  (1) 

Because  of  the  high  H2S04  concentration  (about  5  mol/1) 
in  the  cell  (and  in  the  TLF)  at  100%  SOC,  77ei  =  0.  The  van 
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Fig.  1.  Micrograph  of  an  open  gas  pore  in  NAM  and  a  scheme  of  the  cross-section  through  the  NAM  pore. 


der  Waals’s  component  77vdw  can  be  calculated  by  the  Eq.  (2)  we  obtain  the  following  term  for  the  disjoining 
following  equation:  pressure: 


H  vdw  —  — 


Ah 

6nh? 


(2) 


where  AH  is  the  Hamaker’s  constant  of  the  three  phase 
system  (gas,  electrolyte,  lead).  The  Hamaker’s  constant 
AH  reflects  the  van  der  Waals’s  interaction  between  all 
molecules  of  the  three  phases  in  the  system.  In  our  case, 
Ah  can  be  presented  by  the  following  equation: 


Ah  —  Aei/ei  —  APb/el  (3) 

where  the  Hamaker’s  coefficient  Ael/el  accounts  for  the  total 
interaction  between  the  electrolyte  molecules  only,  and 
Apb/ei  stands  for  the  total  interaction  between  the  electrolyte 
molecules  and  the  lead  crystal  atoms.  It  can  be  approximated 
that  Aei/ei  is  equal  to  the  Hamaker’s  constant  for  water,  which 
is  about  4  x  10“ 13  erg.  APb/ei  (erg)  can  also  be  evaluated 
approximately  using  the  equation 

Apb/ei  «  (Ael/elAPb/Pb)1/2  =  13  x  10-13  (4) 

Apb/pb  is  assumed  to  be  about  40  x  10“ 13  erg.  Hence,  for 
Ah  we  obtain  the  value  of  (—9)  x  10“ 13  erg.  It  must  be 
pointed  out  that  this  value  has  a  negative  sign.  According  to 


Fig.  2.  Schematic  representation  of  the  balance  of  forces  determining  the 
stability  of  the  wetting  film  formed  in  an  open  gas  pore  in  NAM. 


3  x  10"13 

nvdW“  2 nh*  (5) 

The  positive  sign  of  77  means  that  there  is  a  repulsive  force 
between  the  two  surfaces  of  the  TLF,  and  hence  the  film  can 
not  rupture.  Thus,  the  latter  is  thermodynamically  stable. 
Once  formed,  the  wetting  film  will  get  thinner  with  time 
under  the  action  of  the  capillary  pressure  PG  =  2g/R ,  (R  is 
the  radius  of  the  pore  on  whose  wall  the  film  is  formed,  and  a 
is  the  surface  tension  of  the  liquid  phase). 

Ignoring  the  influence  of  the  water  obtained  on  the  film 
thickness,  and  assuming  that  the  film  is  in  a  quasi-equili¬ 
brium  state,  Tlydw  =  Pa  =  2g/R.  So  for  the  thickness  of  the 
TLF  we  obtain: 


,  3/3  x  10“13/? 

*=V  4„ 


<6> 


This  equation  shows  that  h  depends  on  the  open  pore  radius 
R  and  on  the  surface  tension  o .  The  latter  depends  on  the 
presence  of  expanders.  The  thickness  of  the  film  decreases 
with  decrease  of  R.  When  the  film  thickness  becomes  less 
than  twice  smaller  than  the  diameter  of  one  expander 
molecule,  an  additional  repulsive  force  between  the  surfaces 
of  the  film  will  appear.  This  force  is  a  result  of  the  steric 
surface  interactions  between  the  expander  molecules 
adsorbed  at  both  interfaces.  It  stabilizes  the  thickness  of 
the  TLF  [15-19]. 

For  film  thickness  above  100  nm  the  interactive  force 
between  the  two  TLF  surfaces  weakens  rapidly  on  increas¬ 
ing  h  and  the  film  turns  into  a  thick  liquid  layer.  In  the  largest 
open  pores,  the  diffusion  transport  of  02  is  very  slow 
because  of  the  great  film  thickness  and  small  oxygen  con¬ 
centration  gradient  in  it.  Fig.  3  shows  the  dependence  of  the 
differential  pore  volume  and  of  the  pore  surface  as  a  function 
of  pore  radius  for  two  negative  plates  used  in  this  study.  It  is 
evident  that  the  pores  with  the  largest  radii  have  the  smallest 
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02-  gas 
(+  H2  -  gas) 


R/  nm 

Fig.  3.  Mercury  porosimetric  data  for  NAM  sample  from  cell  #2 
(containing  no  additives)  and  from  cell  #5  (with  a  mixture  of  hydrophobic 
and  hydrophilic  additives). 


share  of  the  total  gas-liquid  surface  area.  This  implies  that 
recombination  will  proceed  at  the  lowest  rate  in  those  pores 
due  to  the  larger  TLF  thickness  and  to  the  smaller  pore 
surface.  Since  at  high  electrolyte  saturation  of  the  negative 
plate  only  the  largest  pores  are  open,  the  COC  efficiency  in 
this  case  will  be  very  low. 

On  decreasing  the  plate  saturation,  smaller  and  smaller 
pores  are  opened.  Their  surface  is  large,  the  capillary 
pressure  in  them  is  higher  and  the  TLF  thickness  is  small. 
Consequently,  the  02  concentration  gradient,  and  hence  the 
recombination  rate,  will  increase.  This  is  the  reason  why  at 
lower  electrolyte  saturation  in  the  plate  the  COC  efficiency 
increases  rapidly. 


4.  Processes  taking  place  in  the  thin  liquid  film 

Fig.  4  shows  a  scheme  of  the  process  of  water  production 
at  the  phase  boundary  metal/TLF  according  to  the  electro- 


Fig.  4.  Scheme  of  the  electrochemical  mechanism  of  02  reduction  in  the 
active  centers  at  the  phase  boundary  metal/TLF. 


chemical  mechanism  of  02  reduction.  Oxygen  diffusion 
through  the  wetting  film  on  the  NAM  and  impeded  charge 
transfer  through  the  EDL  are  assumed  to  be  the  rate-limiting 
steps  of  the  oxygen  cycle. 

The  diffusion  of  02  through  the  liquid  film,  the  latter 
considered  as  two  endless  plane  parallel  surfaces  with  a 
distance  h  between  them,  is  described  by  Fick’s  equation: 


dc  d2c 

dt  dx 2 


(7) 


where  D  is  the  diffusion  coefficient  of  02  in  sulfuric  acid 
solution  (in  cm2  s-1)  and  c  is  the  oxygen  concentration  (in 
mol  cm-3)  in  the  TLF. 

Let  us  formulate  the  initial  and  boundary  conditions  of 
oxygen  diffusion.  They  are  schematically  represented  in 
Fig.  5.  We  can  assume  as  an  initial  condition  the  steady 
state  distribution  of  oxygen  concentration  throughout  the 
film,  which  is  a  consequence  of  galvanostatic  polarization 
(/i  =  constant) 


^  ^1  C 2 

c{  0,x)  =  — - x  +  c2 

n 


(8) 


c  i  is  the  oxygen  concentration  at  the  phase  boundary  gas/ 
TLF,  and  c2  is  the  oxygen  concentration  at  the  phase 
boundary  lead  crystal/TLF.  c2  is  not  equal  to  0  due  to  the 


Fig.  5.  Initial  and  boundary  conditions  for  Fick’s  equation.  The  symbol  c 
denotes  the  02  concentration  (in  mol  cm-3)  in  the  electrolyte  film,  x  and  t 
are  the  direction  and  time,  respectively.  The  dashed  line  illustrates  a 
solution  of  Fick’s  equation  at  a  moment  of  time  t  >  0. 
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fact  that  the  kinetics  of  the  oxygen  reduction  process  is 
controlled  by  diffusion  and  impeded  charge  transfer. 

If  the  current  value  increases  from  Ix  at  t  =  0  to  72  at  a 
moment  t  >  0  the  boundary  conditions  are: 

c(t,  0)  =  c4  and  c(t,  h )  =  c3  (9) 

where  c3  is  the  concentration  of  oxygen  at  the  phase 
boundary  gas/TLF,  and  c4  is  the  concentration  of  oxygen 
at  the  phase  boundary  lead  crystal/TLF  at  /2.  According  to 
the  general  Volmer’s  equation  on  increasing  the  current  that 
flows  through  the  cell  the  concentration  c4  should  also 
increase  compared  to  c2.  However,  due  to  the  reduction 
of  oxygen  water  is  produced.  If  the  water  flow  from  the  metal 
surface  is  sufficiently  intense,  the  concentration  change  will 
be  in  the  opposite  direction,  i.e.  c4  will  decrease  compared  to 
c2.  These  boundary  conditions  describe  a  new  steady  state 
distribution  c(x ),  which  will  be  approached  after  a  long 
enough  time. 

In  order  to  solve  the  diffusion  task,  we  use  the  integral 
transformation  of  Laplace  at  time  t.  The  solution  is 
the  Laplace  image  of  the  dependence  of  the  concentration 
on  v  and  s,  where  s  is  the  analog  of  time  in  the  Laplace’s 
space: 


:0,*)  = 


(c4—C2)ch[(x  —  h)  y/s/D]  +  (c3  +  C\)ch{XyJ  sjD ) 


sh(hy/s/D ) 

(C\  -  c2  \ 

\~irx+ci) 


(10) 


The  02  recombination  current  7rec  is  equal  to  the  02  flow 
which  reached  the  phase  boundary  lead  crystal/TLF,  i.e.  to 
the  concentration  gradient  at  this  interface,  which  can  be 
calculated  from  the  derivative  ( dc/dx )  at  x  =  0.  Using  the 
inverse  Laplace  transformation  the  following  equation  is 
obtained: 


free 

4 F 


=  -SD 


C4  ~  C2\ 
y/nDt ) 


(11) 


where  S  is  the  surface  area  through  which  diffusion  pro¬ 
ceeds,  F  the  Faraday’s  constant.  The  equation  is  obtained 
with  the  accuracy  of  an  additive  constant,  since  at  t  — >  00  the 
gradient  (c3  —  c4)  should  be  considered  instead  of  (c\  —  c2). 
(c3  —  c4)  applies  for  the  steady  state  kinetics  at  the  higher 
current.  That  is  why  Eq.  (11)  should  be  represented  as: 


free 

4 F 


=  -SD 


c  3  ~  c4 
h 


-SD 


C4  ~  c  2 
y/nDt 


(12) 


As  is  evident,  the  recombination  current  is  proportional  to 
the  thin  liquid  film  surface  through  which  oxygen  diffuses. 
The  larger  this  surface,  the  higher  the  recombination  current 
and  the  COC  efficiency.  The  second  term  in  this  equation  is 
proportional  to  t~m.  During  the  initial  moments  of  the 
process,  when  t  is  a  very  small  number,  7rec  <  0.  For  this 
case  Eq.  (12)  does  not  have  real  meaning  and  7rec  should  be 
considered  equal  to  0.  After  a  certain  time  z  the  second  term 


of  Eq.  (12)  becomes  smaller  than  the  first  one  and  7rec 
acquires  a  positive  value. 

Eq.  (12)  shows  that  the  recombination  rate  7rec  depends 
largely  on  the  TLF  thickness  h.  The  bigger  the  h  value  the 
smaller  the  oxygen  reduction  rate.  Figs.  1  and  2  show  that  the 
TLF  differs  in  thickness  in  the  different  areas  of  the  pore.  In 
the  areas  where  the  TLF  is  very  thin  (i.e.  h  is  very  small)  the 
recombination  rate  will  be  very  high,  and  in  other  areas 
where  a  thick  TLF  is  formed,  7rec  will  be  very  small  or 
tending  to  0.  This  approach  illustrates  the  fact,  that  oxygen 
reduction  proceeds  actually  in  a  certain  number  of  active 
centers  where  the  TLF  has  a  minimum  thickness.  Thus,  the 
surface  S  should  be  considered  as  an  effective  sum  of  the 
areas  of  the  active  centers  where  oxygen  reduction  takes 
place. 

Let  us  replace  h  from  Eq.  (12)  with  its  equivalent  deter¬ 
mined  for  the  TLF  from  Eq.  (6).  The  following  equation  is 
obtained  for  7rec: 

Vrec  =  =  2.37  x  104SD{c3  -  c4)(/f  -  S(c4  -  c2) 

47y  y  K 

(13) 

where  vrec  means  the  rate  of  02  flow,  which  recombines  to 
H20.  This  is  the  general  equation  for  the  oxygen  reduction 
rate,  which  accounts  not  only  for  the  oxygen  diffusion  rate 
but  also  for  the  properties  of  the  thin  liquid  film. 

For  a  long  enough  time  t  — >  00,  the  second  additive  term 
of  Eq.  (13)  becomes  very  small  and  can  be  ignored.  In  this 
case  the  following  equation  is  obtained  for  the  steady  state 
oxygen  reduction  rate: 

^  =  2.37  x  104SZ)(c3  -  c4)  (14) 

The  value  of  the  steady  state  current  is  determined  by  the 
Pb  crystal  surface  area  on  which  the  reduction  proceeds,  the 
pore  radius  and  the  surface  tension  at  the  boundary  gas/TLF. 
These  parameters  depend  on  the  negative  active  mass 
(NAM)  structure  and  properties,  and  the  kind  of  expander 
used  for  its  preparation.  There  is  also  a  second  group  of 
parameters  characterizing  the  diffusion  process  in  the  thin 
liquid  film.  These  are  the  oxygen  concentrations  on  both 
phase  boundaries  and  the  diffusion  coefficient  of  oxygen  in 
the  TLF. 


5.  Experimental 

Fig.  6  presents  a  scheme  of  the  material  fluxes  in  a  lead- 
acid  cell  with  oxygen  recombination  for  the  case  when  the 
water  decomposition  current  7  is  larger  than  that  of  oxygen 
recombination  7rec  and  part  of  the  oxygen  evolved  on  the 
positive  plate  leaves  the  cell  at  a  rate  7out.  Assuming  that  the 
cell  is  fully  charged,  the  composition  of  the  outgoing  gas 
will  consist  of  02  and  H2  in  a  volume  ratio  1:2  considering 
Faraday’s  law. 
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Fig.  6.  Scheme  of  the  material  transport  processes  involved  in  the  oxygen 
cycle. 


Hence, 


I  =  Ire 


(15) 


when  the  cell  is  charged  galvanostatically  (/  =  constant) 
The  rate  of  the  outgoing  gas  vout  =  /0ut/4E  (in  mol  s-1) 
can  be  estimated  by  combining  Eqs.  (15)  and  (13): 


Vout  =*7= -2.37  x  104SD(c3  -  +  S(c4  -c2)\l  — 


4  F 


R 


nt 


(16) 


As  is  evident  from  Eq.  (16),  vout  is  a  monotonously 
declining  function  of  time.  Two  components  of  the  rate 
vout  can  be  distinguished:  first  a  non-steady  state  one, 
corresponding  to  the  short  period  of  time  immediately  after 
the  current  value  changes  from  A  to  I2\  and  a  second  steady 
state  one,  corresponding  to  the  case  when  t  — >  oo. 


5.1.  Gassing  rate  monitoring  system 


corresponding  to  a  gas  volume  of  100  cm3  at  an  average  rate 
of  2-3  bubbles  per  second.  It  has  been  estimated,  that  the 
average  volume  of  one  bubble  is  0.2  cm3.  The  relative  error 
of  the  method  is  less  than  10%.  The  rate  of  the  outgoing  gas 
v0ut  can  be  calculated  according  to  the  equation: 


Vout 


d«o2  _  KdN 

d  t  d  t 


(17) 


where  no2  is  the  number  of  oxygen  moles  leaving  the  cell 
and  N  the  number  of  the  bubbles  detected  by  the  monitoring 
system.  K  the  constant  of  proportionality  the  value  of  which 
can  be  determined  using  the  equation  for  the  ideal  gas 
PV  =  nRT.  We  assume  that  in  our  case  the  oxygen  partial 
pressure  accounts  for  one  third  of  the  outgoing  gas  pressure 
(which  equals  atmospheric),  i.e.  Po2  =  Patm/3.  In  our 
experiment  the  volume  of  the  outgoing  gas  is  split  into 
equal  bubble  volume  portions,  i.e.  V  =  iWbubbie*  where 
^bubble  is  the  average  volume  of  a  bubble  as  determined 
by  the  gassing  rate  monitoring  system  calibration.  Hence, 
according  to  the  ideal  gas  approach: 

n  rr  APatm  ^bubble  / 1  o  \ 

n0lRT  = - - -  (18) 


Through  differentiating  this  equation  by  time  and 
regrouping  the  parameters  we  obtain 


d N  _  3 RT 

df  UmbbleP  atm 


dn< 


o2 


d  t 


Using  Eq.  (16)  we  come  to: 
3  RT 


d N  _ 

df  UbubbleT5 atm  \  4E 


(  -  2.37  x  104S£>(c3 


(19) 


\  3  /  au 

c4  )\  -B- 


+  Vc4  —  c2)\/  — ; 

Tit  j 


(20) 


In  our  laboratory  a  method  was  elaborated  for  determin¬ 
ing  the  value  of  vout  through  monitoring  the  number  of  gas 
bubbles  leaving  the  cell  per  second  [20].  A  scheme  of  the 
gassing  rate  monitoring  sensor  is  shown  in  Fig.  7.  The 
outgoing  gas  is  lead  through  the  sensor  forming  gas  bubbles 
in  a  tube  filled  with  water.  The  gas  bubbles  are  counted  by  a 
light  source/photo  detector  couple.  When  a  bubble  is 
evolved  it  interrupts  the  light  beam  towards  the  photo 
detector.  The  pulse  is  detected  and  recorded  by  the  PC. 
The  system  is  calibrated  by  counting  the  number  of  bubbles 


The  method  errors  were  minimized  by  integrating  the 
gassing  rate  since  the  moment  of  applying  the  higher  current 
I2.  Thus  the  dependence  of  the  bubble  number  on  t  can  be 
calculated  from  Eq.  (20): 


m  = 


3  RT 


^bubble -P  atm 


I 

4 F 

+ 


3  RT 


^bubble  T^atm 


Replacing: 

3  RT 


a\  = 


a2  = 


and 


I'bubble^atm  ^  2  V  71 


3  RT 


(22a) 


.  ^  -  2.37  x  10 4SD(c3  -  cA{ 7 — 

Vbubble^atm  W  V  'V  R 


(22b) 


t  =  r 


Fig.  7.  Scheme  of  the  gassing  rate  monitoring  sensor. 
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the  cumulative  number  of  bubbles  up  to  the  moment  of  time  t 
is: 

N(z)  =  a2z2  +  aiz  (23) 

Presenting  the  experimental  data  as  N  versus  tl/ 2  (in  sl/2), 
a  second  order  polynomial  regression  can  be  used  for 
experimental  evaluation  of  the  constants  ax  and  a2. 

5.2.  Negative  plates  and  cells 

In  order  to  estimate  the  influence  of  the  surface  tension  at 
the  gas/TLF  boundary  on  the  oxygen  cycle  efficiency  we 
produced  negative  plates  with  two  kinds  of  additives  and  with 
no  additives.  The  plate  dimensions  were  53/65/2.5  mm.  The 
plates  were  tank  formed  and  dry  charged.  Experimental  cells 
were  assembled  with  one  negative  and  two  positive  plates. 
The  separator  was  AGM  (Hollingsworth  &  Vose)  440  g  m2. 
In  order  to  obtain  various  AGM  pore  radii  different  compres¬ 
sions  were  applied  to  two  of  the  cells.  All  cells  were  coupled 
to  the  gassing  rate  monitoring  system.  The  potential  of  the 
negative  plates  versus  a  Hg/Hg2S04  reference  electrode 
was  monitored  in  all  cells.  The  cells  were  filled  to  the  top 
with  1.27  specific  gravity  H2S04  and  left  to  soak  for  30  min. 
After  that  the  excess  electrolyte  was  removed  through  10  min 
draining.  The  cells  were  polarized  on  a  BITRODE  CN  25-12 
testing  module  and,  prior  to  the  COC  tests,  the  cells  were 
conditioned  by  charging  at  0.25  A  for  2  h  followed  by  10 
charge/discharge  cycles  (discharge:  =  0.25  A  to  Ucut-0ff  = 
1.8  VPC;  charge:  7Ch  =  0.25  A,  U\\m  =  2.53  VPC),  charge 
factor  r]  =  118%,  and  another  10  charge/discharge  cycles 
with  the  same  profile  of  charge/discharge. 

6.  Results  and  discussion 

6.1.  Polarization  profiles 

As  established  earlier  using  the  same  technique  [20]  when 
a  current  f  passes  through  a  charged  cell  the  cumulative  gas 
bubble  number  leaving  the  cells  features  an  abrupt  initial 
increase  (A)  and  then  the  rate  of  gassing  decreases  approach¬ 
ing  an  almost  steady  state  value,  as  shown  in  Fig.  8.  This 
initial  rise  was  attributed  to  formation  of  gas  channels 
through  the  AGM  separator. 

In  this  investigation  our  attention  is  focused  on  the 
phenomena  that  proceed  in  the  negative  plate  after  the 
gas  channels  are  formed  in  the  separator.  For  the  purpose 
of  this  study,  the  cells  were  first  polarized  galvanostatically 
until  a  steady  state  gassing  rate  was  obtained,  i.e.  all  gas 
channels  in  the  AGM  separator  were  open.  Then  a  new 
constant  current  value  was  applied  (period  (B)  in  Fig.  8),  and 
the  changes  in  gassing  rate  were  recorded.  We  assume  that 
the  application  of  the  new  current  does  not  cause  any 
significant  changes  in  the  gas  channels  structure  and  the 
only  phenomena  that  are  affected  by  the  current  change  are 
those  that  take  place  in  the  negative  active  mass  (NAM), 


Fig.  8.  Cumulative  number  of  bubbles  vs.  time  at  polarization  currents  /, 
and  I2  (I i  <  h)- 

being  related  to  oxygen  reduction.  These  can  be  estimated 
indirectly  from  the  changed  rate  of  the  gas  flow  leaving  the 
cell.  This  change  depends  on  the  degree  of  saturation  of 
NAM  with  electrolyte  and  on  the  new  current  value  /2. 

In  order  to  verify  experimentally  the  validity  of  Eq.  (20) 
we  used  two  saturation  levels  and  two  /2  current  values. 

H.  High  electrolyte  saturation  level  of  the  plates  and  of  the 
AGM  separator,  and  a  substantial  current  rise  from 
h  =  50  mA  to  I2  =  150  mA.  In  this  case,  the  increased 
02  pressure  in  the  separator  at  current  12  forces  the  gas  flow 
to  penetrate  into  new  pores  of  the  negative  plate.  They  have 
smaller  radius  R  and  when  filled  with  oxygen  they  take  part 
in  the  recombination  process.  Hence,  the  reaction  surface  S 
increases. 

L.  Low  electrolyte  saturation  level  of  the  active  block  and 
but  a  small  current  rise  from  I\  =  450  mA  to  I2  =  500  mA. 
In  this  case,  no  new  gas  pores  are  opened  in  NAM  and  R  and 
S  remain  unchanged. 

The  investigation  was  performed  using  four  kinds  of 
negative  plates  with  various  surface  tensions  (er)  at  the 
gas/electrolyte  boundary.  The  following  NAM  additives 
were  used: 

(a)  hydrophobic  polymer — cell  #6 

(b)  hydrophilic  polymer — cell  #4 

(c)  a  mixture  of  hydrophobic  and  hydrophilic  polymers — 
cell  #5 

(d)  no  NAM  additives  but  different  AGM  compressions: 
20%  in  cell  #1  and  25%  in  cell  #2 

The  cells  were  assembled  into  a  battery  and  the  latter  was 
set  to  overcharge  and  COC  studies  were  performed. 

6.2.  COC  in  cells  with  high  H2S04  saturation  of  the 
negative  plate 

Let  us  assume  that  Vbubbie  =  constant  irrespective  of  the 
gassing  rate,  Dq2  =  constant  since  the  cell  is  fully  charged 
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and  T  =  constant.  In  this  case  for  Eq.  (20)  we  obtain: 


d N  I  ,  N  „  ,* 

2(C4_C2)~a/t 

-K3(C3-C4)  —  ?- 
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When  the  current  value  is  changed  from  I\  =  50  mA  to 
I2  =  150  mA  the  parameters  R ,  S,  (c3  —  c4)  and,  though  less 
so,  (c4  —  c2)  will  also  change,  o  is  constant  for  this  particular 
cell. 

Fig.  9  presents  the  number  of  bubbles  leaving  the  cell  as  a 
function  of  t  for  all  five  cells.  What  conclusions  can  be 
drawn  based  on  the  data  in  Fig.  9  and  Eq.  (24)7  All  curves 
rise  initially  and  approach  a  steady  state  value.  This  indicates 
that  after  a  certain  non- steady  state  period  of  time  the  COC 
reaches  a  certain  efficiency  and  gets  into  a  steady  state  stage. 
What  are  the  changes  during  the  non-steady  state  period 
according  to  Eq.  (21)7  On  raising  the  current  the  oxygen 
pressure  at  the  boundary  AGM/NAM  increases.  The  gas 
penetrates  into  the  pores  with  smaller  radii  (replacing  the 
electrolyte)  and  thus  the  reaction  surface  S  increases. 
According  to  Eq.  (13),  the  increased  S  surface  accelerates 
the  rate  of  02  reduction.  The  02  concentration  gradient 
(C3  —  c4)  between  the  two  TLF  interfaces  increases.  Hence, 
the  third  term  in  Eq.  (24)  will  acquire  a  value  close  to  those 
of  the  other  two  terms.  The  difference  between  them  deter¬ 
mines  the  gassing  rate  dN/dt.  The  derivative  decreases 


Fig.  9.  Dependence  of  cumulative  number  of  bubbles  on  time  at 
I  =  150  mA. 


Fig.  10.  Dependence  of  cumulative  number  of  bubbles  on  the  square  root 
of  time  at  I  =  150  mA  (solid  lines)  and  results  from  the  fit  by  Eq.  (21) 
(symbols). 

approaching  0.  This  is  exactly  the  form  of  the  experimental 
curves  in  Fig.  9. 

What  is  the  duration  of  the  non- steady  state  stage  and 
when  will  the  steady  state  stage  begin?  This  can  be  estimated 
by  fitting  the  experimental  data  to  Eq.  (23). 

The  number  of  bubbles  leaving  each  cell  is  shown  in 
Fig.  10  as  a  function  of  the  square  root  of  time  (s).  The 
experimental  data  were  fitted  by  a  second  order  polynomial 
regression,  which  corresponds  to  Eq.  (23),  where  N  is  the 
number  of  bubbles  leaving  the  cell  for  the  period  of  time 
between  0  and  t.  The  coefficients  ax  and  a2  are  related  to  the 
parameters  presented  in  Eq.  (22).  The  experimental  data  are 
illustrated  by  a  solid  line  in  the  figure,  and  the  best  fit  values 
calculated  according  to  Eq.  (21),  by  symbols.  The  symbols 
used  for  the  first  and  the  second  stages  are  different,  and  the 
transition  point  is  marked  by  a  vertical  line.  The  equation  fits 
well  all  parts  of  the  experimental  data  plot,  the  coefficients 
a i  and  a2  have  different  values  for  the  two  stages.  Their 
values  for  the  non-steady  state  (I)  and  for  the  steady  state  (II) 
stages  are  given  in  Table  1. 

According  to  Eq.  (22a)  the  increase  of  a x  in  the  second 
stage  is  due  to  the  increase  of  the  reaction  area  S  where 
oxygen  reduction  proceeds.  The  data  in  Fig.  3  show  that  this 
area  increases  rapidly  in  the  pores  with  smaller  radii.  The 
oxygen  concentration  difference  (c4  —  c2)  at  the  boundary 
metal/TLF  also  increases. 

According  to  Eq.  (22b),  the  fact  that  a2  has  a  negative  sign 
during  the  second  stage  indicates  that  S  and  (C3  —  c4) 
increase  and  R  decreases.  Hence  2.37  x  10~4SD(c3  —  c4)~ 
(ncr/R)1/3  >  (//4T1)  and  a2  <  0. 

Fig.  9  shows  that  cell  #6  with  hydrophobic  and  cell  #4 
with  hydrophilic  polymer  additives  have  the  lowest  COC 
efficiency.  For  cell  #6,  the  transport  of  H+  ions  in  the  pores 
of  NAM  is  impeded  by  the  hydrophobic  layer.  In  the 
negative  plate  of  cell  #4  the  gas  pores  are  blocked  by  the 
hydrophilic  polymer  and  02  reaches  only  the  largest  pores 
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Table  1 

ai  and  a2  values  for  the  polynomial  regression  in  the  case  of  high  H2S04  saturation  level  of  the  negative  plate 


Current 

Stage 

Coefficient 

Cell  1 

Cell  2 

Cell  4 

Cell  5 

Cell  6 

0.15  A  (0.09  C) 

1st  stage 

ai  [s“1/2] 

2.97 

2.37 

8.05 

0.98 

10.19 

a2  [s-1] 

0.070 

0.073 

0.071 

0.092 

0.052 

2nd  stage 

ai  [s'172] 

28.06 

9.56 

41.48 

20.54 

116.96 

a2  [s-1] 

-0.149 

-0.053 

-0.179 

-0.107 

-0.520 

which  have  but  a  small  surface.  When  both  polymer  addi¬ 
tives  are  used  in  combination  the  obtained  balance  between 
the  hydrophobic  and  hydrophilic  properties  of  the  surface 
causes  a  dramatic  improvement  in  COC  efficiency,  which  is 
observed  in  cell  #5. 

The  cells  with  no  additives  (cell  #1  and  cell  #2)  have  high 
COC  efficiency.  The  latter  is  higher  for  the  cell  with  higher 
compression,  since  the  02  pressure  inside  the  NAM  pores  is 
larger.  The  gas  penetrates  into  the  pores  with  smaller  radii 
and  the  reaction  surface  S  increases.  Consequently,  the  COC 
efficiency  improves.  An  interesting  fact  from  a  practical 
point  of  view  is  that  the  cell  with  a  mixture  of  hydrophobic 
and  hydrophilic  polymer  additives,  and  20%  compression 
(cell  #5)  has  a  higher  COC  efficiency  than  the  cell  under  the 
same  compression  but  with  no  additives  (cell  #1).  This 
indicates  that  a  proper  selection  of  NAM  additives  with 
hydrophobic  and  hydrophilic  properties  can  influence  both  a 
and  h  (film  thickness),  and  hence  the  COC  efficiency. 

6.3.  COC  in  cells  with  low  H2S04  saturation  level  of 
the  negative  plate 

The  current  value  was  changed  from  450  to  500  mA.  In 
this  case,  the  pores  with  small  radii  are  involved  in  the 
oxygen  reduction  process  during  the  initial  stage  already 
(450  mA).  The  current  rise  does  not  affect  the  values  of  R 
and  S.  Since  the  current  change  is  relatively  small  the 
pressure  does  not  change  either,  a  remains  constant,  too. 
In  this  case  Eq.  (21)  is  reduced  to: 

^  =  K[I  +  K'2(c4  -  C2)  Y7  -  K'3(c3  -  c4) yz  (26) 

Fig.  1 1  presents  the  experimental  data  obtained  as  number 
of  bubbles  versus  polarization  time  for  l  —  500  mA.  When 
the  saturation  is  low,  the  steady  state  stage  is  reached  very 
quickly  and  the  figure  illustrates  mainly  this  stage. 

An  analysis  of  the  experimental  curves  considering 
Eq.  (26)  leads  to  the  following  conclusions.  When  the 
current  rises  from  l\  to  /2,  the  02  concentration  difference 
at  the  two  TLF  surfaces  increases  due  to  the  02  concentra¬ 
tion  rise  at  the  boundary  pore/TLF.  This  causes  the  oxygen 
reduction  rate  to  increase.  Water  is  formed  and  the  TLF 
thickness  grows  causing  a  decrease  of  the  rate  of  oxygen 
diffusion  towards  the  metal  surface  and  hence  a  decrease  of 
the  electrochemical  reaction  rate.  Under  the  action  of  the 
capillary  forces  the  TLF  tends  to  get  thinner.  So  after  a 


certain  period  of  time,  a  dynamic  stage  is  reached  when  the 
differences  (C3  —  cf)  and  (q  —  c2)  become  constants.  On  the 
other  hand,  the  second  term  in  Eq.  (26)  becomes  very  small 
with  time.  It  can  be  ignored  because  its  value  is  much 
smaller  than  that  of  the  first  term.  So  dN/dt  =  constant. 
This  is  actually  observed  in  Fig.  11,  i.e.  a  linear  increase  of 
the  cumulative  number  of  bubbles  with  time. 

Fig.  12  gives  the  dependence  of  the  number  of  bubbles 
versus  tl/2.  The  experimental  data  are  illustrated  by  a  solid 
line,  while  the  best  fit  values,  calculated  according  to 
Eq.  (21),  are  denoted  by  symbols.  The  experimental  data 
fit  well  to  Eq.  (21)  within  the  entire  time  range  (12  h).  The 
low  electrolyte  saturation  of  the  NAM  yields  a  system  of  gas 
pores  with  a  large  surface  and  a  thin  TLF.  This  ensures  high 
COC  efficiency  for  a  long  period  of  time. 

Table  2  presents  the  values  for  ax  and  a2  for  the  best  fits 
calculated  according  to  Eqs.  (21)  and  (23). 

Though  the  surface  S  is  substantially  larger  in  this  case, 
the  ax  values  are  lower  than  those  in  Table  1  (Eq.  (22a)). 
Probably,  when  the  higher  current  /2  is  applied,  the  02 
concentration  on  the  lead  surface  does  not  change  substan¬ 
tially  and  the  concentration  difference  (q  —  c2)  is  rather 
small.  The  value  of  ax  remains  very  small. 

The  greater  current  value  I  =  500  mA  makes  the  first  term 
in  Eq.  (22b)  greater  than  the  second  one  and  hence  the  sign 
of  a2  is  positive. 


Fig.  11.  Dependence  of  cumulative  number  of  bubbles  on  time  at 
/  =  500  mA. 
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Fig.  12.  Dependence  of  cumulative  number  of  bubbles  on  the  square  root 
of  time  at  I  =  500  mA  (solid  lines)  and  results  from  the  fit  by  Eq.  (21) 
(symbols). 


Table  2 

a i  and  a2  values  for  the  polynomial  regression  in  the  case  of  a  low  H2S04 
saturation  level  of  the  negative  plate 


Current 

Coefficient 

Cell  2 

Cell  4 

Cell  5 

Cell  6 

0.50  A  (0.31  C) 

ai  [s'172] 
a2  [s-1] 

1.29 

0.0112 

1.17 

0.0121 

1.63 

0.0135 

1.51 

0.0220 

The  experimental  data  discussed  above  support  the  con¬ 
clusions  drawn  from  Eq.  (21)  thus  providing  an  experimental 
proof  of  the  model  for  oxygen  diffusion  and  impeded  charge 
transfer  through  a  thin  liquid  film  as  the  limiting  stages  of  the 
processes  involved  in  the  closed  oxygen  cycle. 

7.  Conclusions 

The  aim  of  the  present  work  was  to  study  the  processes 
taking  place  in  the  thin  liquid  film  (TLF)  that  covers  the  Pb 
crystals  of  NAM.  Oxygen  diffusion  through  the  TLF  cover¬ 
ing  the  NAM  metal  surface  and  impeded  charge  transfer 
through  the  electric  double  layer  were  assumed  to  be  the 
rate-limiting  steps  of  the  oxygen  cycle.  These  processes 
depend  on  the  TLF  properties,  the  structure  of  NAM  and  on 
the  electrolyte  saturation  level  of  the  plate. 

Based  on  the  colloid  systems  theory  it  has  been  estab¬ 
lished  that  the  thickness  of  the  TLF  depends  on  pore  radius 
and  on  the  surface  tension  at  the  gas/liquid  interface. 


Solving  Fick’s  equation  for  the  defined  initial  and  bound¬ 
ary  conditions  we  obtained  an  equation  giving  the  recombi¬ 
nation  current  as  a  function  of:  the  metal  surface  area  on 
which  recombination  takes  place,  oxygen  diffusion  coeffi¬ 
cient,  thickness  of  TLF,  oxygen  concentration  at  the  two 
TLF  surfaces  and  polarization  time. 

Substituting  the  value  for  the  TLF  thickness  by  an  expres¬ 
sion  based  on  colloid  chemistry  theory,  a  general  equation 
was  obtained  for  the  relation  between  oxygen  recombination 
rate  and  the  properties  of  the  TLF,  NAM  structure  and  time. 
The  validity  of  this  equation  was  proved  experimentally  for 
low  and  high  electrolyte  saturation  levels  of  the  plates  and 
for  negative  plates  containing  hydrophilic  and  hydrophobic 
additives  to  NAM,  as  well  as  for  two  active  block  compres¬ 
sions.  The  experimental  results  are  in  perfect  agreement  with 
the  conclusions  drawn  based  on  the  model  for  the  slowest 
stages  of  COC  proposed  in  this  work.  This  makes  the  model 
applicable  to  drawing  practical  conclusions  about  the  meth¬ 
ods  to  be  used  for  improving  the  efficiency  of  the  COC,  and 
about  the  future  trends  in  developing  advanced  negative 
plate  technologies  for  VRLA  batteries. 
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